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Abstract 

We consider the Riemann problem for relativistic flows of polytropic fluids and find relations for 
the flow characteristics. Evolution of physical quantities take especially simple form for the case of 
cold magnetized plasmas. We find exact, explicit analytical solutions for one dimensional expansion 
of magnetized plasma into vacuum, valid for arbitrary magnetization. We also consider expansion 
into cold unmagnetized external medium both for stationary initial conditions and for initially 
moving plasma, as well as reflection of rarefaction wave from a wall. We also find self-similar 
structure of three-dimensional magnetized outflows into vacuum, valid close to the plasma-vacuum 
interface. 

The key results of this work, the self-similar solutions, were incorporated post-initial submission 
into appendices of the published version of Granot et al. (2010). 



I. INTRODUCTION 



Relativistic shock waves are common in different physical systems [T], from heavy ion 
nuclear collision [e.g. |^ to astrophysical shocks in pulsar winds [3j, Active Galactic Nuclei 
[e.g. E] and Gamma Ray Bursts [S]. Many modern computational algorithms are based on 
the solution of Riemann problems [e.g. 6J. These algorithms are based on Godunov-type 
shock-capturing schemes and do not require large artificial viscosity or smoothing operators. 
Analytical solutions to the corresponding Riemann problems are then important for code 
testing. 

Exact, explicit non-linear solutions of relativistic fluid equations, and especially relativis- 
tic MHD equations, are rare. In a general form the relativistic Riemann problem was solved 
by O [8] , who find the solutions for Riemann invariants and for the characteristics in quadra- 
tures. In this paper we find simple expressions for the characteristics of Riemann simple 
waves, and in particular for the astrophysical important case of dynamics of a cold, relativis- 
tically magnetized plasma. Our results can be used for benchmark estimates of the overall 
dynamical behavior in numerical simulations of relativistic flows and strongly magnetized 
outflows in particular. 

II. ONE DIMENSIONAL EXPANSION OF POLYTROPIC GAS INTO VACUUM 
A. Polytropic equations of state 

Let us assume that pressure is a polytropic function of density, P = Kp^, with constant 
F; then the internal energy density, excluding rest mass is S = P/(r — 1), enthalpy h = 
p + £ + P = p + T / {T — 1)P and sound speed can be expressed as 




pdP (T-l)TP/p 



(1) 



hdp (r - 1) + TP/p 



The corresponding four-velocity is 




(r-i) + r(2-r)p/p 



(r - i)TP/p 



(2) 
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B. Riemann invariants 



Consider a one-dimensional flow of fluid along z direction, neglecting variations of quan- 
tities across the flow. The governing equations are 

5t(7P)+5.(7/3p) = 
S^Too + 9,To, = 
9iTo, + d,T,, = 
Too = l\£ + P + p)-P 
To. = 'y'P{S + P + p) 

T,. = l''(3\8 + P + p) + P (3) 

where P is pressure, £ is energy density (excluding rest mass), p is density and (3 and 7 are 
fluid's velocities and Lorentz factors. 

Assume that initially the fluid occupies region z < with constant density and pressure 
and that expansion proceeds into positive z direction. Fig [T} 

Combining the conservation of energy and momentum, we find [9l [10] 

{dt + 13d,) (3 = - , p^ 2 

[t + p + Pj'y^ 

{dt + (3d;) {S + p) = -{S + p + P)7' {(3dt + d-,) (3 (4) 

These two equations and the equation of state combine to give equations for Riemann 
invariants 

J.^llnl±4±-^lnfII + -- (5) 

2 1-/3 -cs ^ ' 

[?]. (J+ corresponds to the forward characteristic). The combinations {(3 ± /3cJ/(l ± (3c^(3) 
can be recognized as relativistic transformation of velocities, plasma velocity and wave phase 
velocity. Riemann invariants are constant along characteristics 

13 ±13,. 

= r±M: '''' 

(the forward characteristics correspond to upper sign). 
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FIG. 1. Geometry of the flow. Initially the fluid at rest occupies region z < with constant density 
and pressure and that expansion proceeds into positive z direction. The fluid moves in the positive 
z-direction while the front of the rarefaction wave propagates in the negative z direction. 

Equations ([s]) simplify if we introduce Doppler factor 

(Under co-linear Lorentz transformation, the corresponding Doppler factors change as 5 ^ 
(J'^boost) where 6' is a Doppler factor in a frame moving with respect to the observer frame 
with Doppler factor 5boost- Doppler factors multiply for aligned velocities or divide in case 
of counter aligned velocities.) 

Next we assume that all quantities depend on self-similar coordinate rj = z/t. When 
expressed in terms of Doppler factors, Eqns. (|5| then simplify 

{SlSl - 5l)d,J+ = (8) 
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(7) 



This gives 5^ = 5j^5c^ for the forward and 5p = 5r;/^cs for the backward characteristics 
corresp ondingly. 

Eq. (|9| is identically satisfied on the forward characteristics, while Eq. ^ is satisfied on 
the backward characteristics. Using the expressions for Riemann invariants and the shape 
of characteristics we find 

d,iJ.i5^ = 5J5J) = Q (10) 

Which gives the full solution for simple waves in relativistic fluids: 

C / 1 — c / /T — 1 \ ^^^^ 
Sri + = -r^ \ — ) , = 5c Sr., forward characteristics (11) 

dc, \i + c,/vr - 1/ 

(I -\- c j i/r — 1 \ ^"^^^ 
^ . , (5/3 = 5v/8csi backward characteristics (12) 
1 - c^/vr - 1/ 



where Cq is a constant to be determined from initial condition. Relations ( 12 ) give a general 
relation between the self-similar coordinate 77, local sound speed in the flow Cg and velocity 
of the flow /3. These are transcendental equations for Cs (77) and P{ri) 



C. Magnetized cold plasma 

In many astrophysical phenomena the magnetic field controls the overall dynamics of 
plasma. These likely include magnetars (strongly magnetized neutron stars possessing super- 
strong magnetic fields), pulsars and pulsar wind nebulae, jets of Active Galactic Nuclei and 
Gamma-Ray Bursters. The plasma of these exotic objects can be described as relativistically 
strongly magnetized. This means that the inertia of this plasma is dominated by the mag- 
netic field and not by the particle rest mass, -B^/Svr ^ pc^, and that the propagation speed 
of Alfven waves approaches the speed of light. Thus, the conditions in such plasmas are 
very different from the conditions encountered in laboratory plasmas, plasmas of planetary 
magnetospheres, and the interplanetary plasma. 

For example, magnetic fields may play an important dynamical role in the Gamma 
Ray Burst outflows [e.g. [HI [12]. They may power the relativistic outflow through, e.g.., 
Blandford-Znajek [13] process [see also 111] . and contribute to particle acceleration in the 
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emission regions. The extreme physical conditions of magnetically dominated plasma suggest 
new physics and demand a systematic study. In relativistic magneto-plasma Godunov-type 
schemes are discussed by [T5] . 

Consider cold ideal magnetized plasma moving perpendicular to the direction of the 
magnetic field. This case, in fact, is just a particular case of polytropic equations of state 
P (X for unmagnetized fluid flow [16]. The case of cold magnetized plasma is, in fact, 
somewhat special since it gives simple relations between density and sound four-velocity, as 
we demonstrate below. For cold plasma, when the pressure is purely magnetic, P = 
(for convenience, we renormalize magnetic field by ^/4^l below), polytropic index F = 2, and 
defining magnetization parameter a = B"^ / p, Eq. ^ gives 

Us = Ua = -j^= = (13) 

the four- velocity of Alfvenand fast magnetosonic waves in strongly magnetized plasma [31 
(3a is Alfven velocity] . 

In case of magnetized plasma, the equations of motion ^ should be supplemented with 
induction equation and definitions of magnetic energy density and pressure 

9t(75) + 9,(7/35) =0 

£ = P = (14) 



where B is a proper magnetic field. From Eqns (14), it immediately follows that magnetic 
field and density are proportional, B = {p/po)Bo, where po and Bq are constants taken to 
be density and magnetic field in the initial state. 
For F = 2, Riemann invariants take the form 



J+ = log (S^Sl) 

J_ = log {5^5 f) (15) 



where we introduced Doppler factors Sa 



where (3a is the local Alfven speed. 



General solution (12) then gives 

6a = Co5f /^ 6p = 6r,6f = Co5l'^ (17) 
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where Cq is a constant to be determined from initial condition. Upper signs in Eq. (17) 
correspond to the forwards characteristics. 



1. Stationary initial conditions 

Assuming that initially the plasma is at rest, so that at the front of the rarefaction 
wave (where /3 = 0, 5^ = 1) the Alfven velocity corresponds to the Alfven velocity in the 
unperturbed medium, given by corresponding Doppler factor 6 Aft, we find a simple fully 
analytical solution to the problem of one dimensional relativistic expansion of magnetized 
gas into vacuum, see Fig. |2| 
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2/3 



= ^5'/^, 5^ = ^ (18) 

Or, 

The other solution, corresponding to expansion to the left is 

A2/3 

= ^,,^A = 6'Il6l/\ (19) 

These solutions give the velocity and magnetic field as functions of the self-similar coordinate 
T] and initial magnetization a. 

The relations for the more commonly used parameters are (for the forward characteristics) 



7 = o '^;'^i,o 



e2/3 r2/3 ^ 

"4 n -r 2/3 r2/3 



1 + sTX" 



u 



A 




B 



I^A- ^ - 2/3 .4/3 
^'7 + '^A.O 



I + Paa 

U,0 ' 



^^_^^„,/3a,o = V^/(l + ^) (20) 



Plasma density at each point is p = {U\/a)pQ and magnetic field is -B = {U\/a)Bo. We 
stress that these solutions are exact, no assumptions about the value of the parameter a 
were made. 
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In the strongly magnetized limit o" — )■ oo we find 

1 



/3 = 1 



7 



21/354/3^2/3 

2\ 1/3 



a 



1/3 



cr 

^vj V 27 

The front of the rarefaction wave is located where dp = 1. This gives 6rj = 1/6a,o, 



26„ 



a 



Vrw 



1+CT 



(21) 



(22) 



Thus the rarefaction wave propagates into undisturbed plasma with Alfven velocity in the 
unperturbed plasma (3a = 



l + (7 
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FIG. 2. Four- velocity of the flow and Alfven four velocity Ua for one-dimensional self-similar 
expansion of magnetized gas into vacuum. For this plot a = 2. 



The vacuum interface corresponds to zero magnetic pressure, where (5a = 1, Eq. (20). 



This occurs at 5^ = 5vac = S\ o (this corresponds to the maximum allowed value of 5,,). 



Vacuum interface is located at 



(l + 2a) ^ ifa«l ^ ^ 



The Lorentz factor of the vacuum interface is exactly 



7,,, = = 1 + 2(T (24) 

'I - T]^ 



/vac 



The corresponding four-velocity is 



/ / , , 2o- ifo- > 1 

t^vac = '7vac7vac = 2V0-(1 + O") ^ <( (25) 

2v/(T ifo- < 1 



Since in the limit a <^ 1 the initial Alfven velocity is Pa,o = \^<^, the front expands with 
fvac = 2/3yi,0; in agreement with known results for polytropic fluid [expansion into vacuum 
with velocity 2c,/(r - 1) [nHH] 

The flow becomes sonic, Pa = P always at rj = 0, 6^^ = 1; this is similar to the general 
case of non-relativistic expansion of polytropic gas into vacuum [201 §20]. At this point the 
four-velocity of the flow 

^ f (a/2)V3 ifoi 
U{v = 0)^{' '' (26) 
[ 2v/^/3 ifa < 1 

Equation for forward and backward characteristics are 



dtz\ 



l + m 1 + 5]ll5T 



The backward characteristics is a straight line, z = rjt, with rjrar < f] < rjvac- Equation for 
the forward characteristics integrates to give 

^ = '^"^' S .8/3^^.4/1 3/2 ' ^'^ < ^Xo (28) 



For a characteristics that starts at a point t = to, z = —toy a /{I + a), this gives 

^^(«,o-l)v/^^^|^ (29) 

Thus, the forward characteristics never cross the vacuum interface (5^ = S\q) and be- 
comes parallel to it as t — t- oo. 



The flow lines are given by 



Which integrates to give 



* = Const^^^^ii*|^,(5, <(5^„ (31) 



In particular, for a flow line starts at a point t — to, z — —toJcT/(l + a), this gives 

t _ {sXo - (1 + s',) 



(32) 



^0 i+sxo {sfo-st^"r 

Thus, both the characteristics and the flow lines asymptote to — )> 6\q. 

For 0" ^ 1 the expansion is relativistic practically in the whole flow. In the bulk, 7 ~ a^^^ 
(e.g., at the sonic point 77 = 0, 7 = (cr/2)^/^). In a narrow region near Tyyao with thickness of 
the order A77 ~ the Lorentz factor of the flow approaches 7vac (at this A77 the Lorentz 

factor reaches half of its maximum values). As the flow expands, the local magnetization, 
UJ^, decreases. 

The energy flux, 

r.. = T-(p + S)/J= ^^^^ p. (33) 

reaches maximum at = 0. In the limit o" ^ 1 it stays nearly constant in between rjnw and 
rjyac at a value Tq^ ~ -Bo/4. The energy density component, 

T,„ ^,(B+p)-B/2 = ^^^^ p„ (34) 

decreases towards the vacuum interface, remaining nearly constant in between rjRw and rjyac 
at the same value Too ~ -Bo/4. It may be verified by direct calculations that the total energy 
in the flow, integral of Too from rjjm/ to ^^^aci equals the total energy in the initial state 
between rjnw and 0, with the energy density equal to po(l + cr/2). 



2. Moving piston: expansion into vacuum 

Let us now assume that in the undisturbed plasma, the piston, is moving with velocity (3^ 
towards the external medium. The corresponding relations can be trivially obtained using 
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Riemann invariants and imposing a condition that the piston is moving with velocity j3w 
and corresponding Doppler factor 5^ — ^^(1 + f3w)/{l — Thus, in the previous relations 
we need to make a substitution Saa — > (^a.oa/^ (this choice assumes that the magnetized 
medium moves toward the unmagnetized one; in the opposite case, 6a,o ~^ Sa,o/V^)- We 
then find exact solution of relativistic Riemann problem for expansion into vacuum of cold 
strongly magnetized plasma with magnetization parameter a and moving with velocity (3^,: 



1 + 5T.Sr6'J' 



^2/3 

(In case of a boost away from the contact, we should substitute — >■ 5~^/^, in accordance 
with the Lorentz transformation of the Doppler factor.) 

The vacuum interface is moving with (5^ — which in the limit cr, 7^,, ^ 1 gives 

7mc = 27^(1 + 2(7). 

The front of the rarefaction wave is located aX 5^ — 5^, which gives 

5, = -^^^ (36) 

The front of the rarefaction wave is stationary (located at 77 = 0, 5^ = 1) when the flow 
moves sonically, with the Alfven velocity of the undisturbed plasma, = /3a,o- We stress 
that these solutions are exact, no assumptions about the value of the parameter a and the 
velocity jiy, were made. 



3. Reflection of the rarefaction wave from the wall 

Let us now assume that the initial state with — Q occupies a limited region of space 
— L < 2; < 0, with impenetrable wall at 2; = — L. As the rarefaction wave reaches the 
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wall, it will be reflected, creating a secondary RW. The initial RW propagates to the left, 
while the secondary RW propagates to the right into the region disturbed by the first RW. 
Let's denote f3Rw,2 the velocity of the secondary RW. Its front propagates along the forward 



characteristics with velocity given by Eq. (29), Fig. 3 



In this case the expansion is non-self-similar, since there is a typical scale in the problem 
L. In addition, even at asymptotically long times after the beginning of the fiow, at times 
t ^ L/c, the expansion of a magnetic shell is not self-similar, since conservation of mass 
(and magnetic fiux) results in different scaling of magnetic (oc B'^) and rest mass energy 
densities. 

Still, it might be possible to find exact solution for the expansion of a magnetized layer into 
plasma following the corresponding hydrodynamic approach [21]. (One can always replace 
the nonlinear problem of one-dimensional hydrodynamic evolution with a linear equation 
using [22] transform). We leave this problem for future; here we consider the motion of the 
front of the secondary RW. 

Let us discuss the limiting cases of fiow evolution long after the refiection of the secondary 
RW from the wall, t ^ Iq. In the high magnetization limit, 5a,o = 2-\/o" ^ 1, and for large 



Lorentz factors of the the secondary RW, I3rw,2 ~ 1 — 1/(27|.^^2)) Eq. (29) gives 

,1/3 

(37) 



2 

RW,2, 



For a long time, when t < 2to7RVK2? the rarefaction accelerates only slowly, 

7RW,2 = V2a^'\t/tof'\ iH < 8toa^ (38) 



For t > 2to7^iy25 which using Eq. (38) implies t > Stocr , the secondary RW reaches 
the terminal Lorentz factor 'yRw,2 ~ 2(j, the Lorentz factor of the vacuum interface. In 
self-similar coordinates, the front of the secondary RW moves from r] = —\J<y/ (1 + o") to 
r] = 2^(7(1 + (t)/(1 + 2(t), Fig. jsj Note, that even in the case when the second rarefaction 
wave approaches the vacuum interface, the expansion is still generically non-self-similar 
since there are two conserved quantities, energy and magnetic fiux, that scale differently 
with radius. 



The velocity of plasma in front of the refiected rarefaction waves is given by Eq. (30). 
The fiow initially accelerates 

7 = {t/t,f\"' (39) 
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and reaches terminal Lorentz factor 7 = 1 + 2cr at times t > Stocr'^ (see also [? ]). As the 
secondary rarefaction wave catches up with the flow, the flow is decelerated. Thus, most 
of the acceleration occurs in a regime when the flow is causally disconnected from the wall 
(contrary to the claim in Ref. [23], version 1). 




FIG. 3. Left Panel: Lorentz factor of the secondary rarefaction wave as function of t' , time 
since reflection. The secondary RW propagates in the flow disturbed by the primary RW. In the 
limit t' ^ Sfj^to) the secondary RW reaches the terminal Lorentz factor of the vacuum interface 
7vac = 1 + 2(7 (dashed Une). The plot is for a = 2. The limiting value of the Lorentz factor is 
reached very slowly, with jiiw,2 = V^cr'^^^it/tQy/^ for t ^ Sa^tQ. Right Panel: motion of the 
secondary RW in self-similar coordinates r] = z/t. At time to = secondary RW is launched 
from z = —\/a/{l + cr)to, corresponding to r] = —y^a/{l + a) (lower dashed line) in the forwards 
direction. At t — )■ 00 the front of the RW approaches the vacuum interface rj = 2^a(\ + a)/ {l + 2a) 
(upper dashed line). 

III. ONE-DIMENSIONAL EXPANSION OF MAGNETIC PISTON INTO PLASMA 
A. Stationary initial conditions 

The solution for the expansion of magnetized plasma into vacuum derived in ^|Tl] can be 
used to consider a decay of a contact discontinuity between magnetized plasma with density 
Po and magnetization parameter a and a cold plasma with density pex- For clarity, we refer 
below to the magnetized component as a piston, and the non-magnetized plasma as the 
external medium. 
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For cold external plasma, we expect the formation of four regions: undisturbed external 
plasma, shocked external plasma, expanding magnetic piston and undisturbed magnetized 
plasma. These regions are separated by the forward shock, a contact discontinuity and an 
expansion wave. We expect that the dynamics of the system will be self-similar. 

The solution for simple waves derived in ^|TT] can be used to describe the expansion wave 
part of the flow. For non-zero external density, the simple waves will terminate at a contact 
discontinuity (CD), located at a particular value of the self-similar variable rjcD < Vvac, 
so that the magnetic pressure of the expansion wave balances the thermal pressure of the 
shocked medium, which, in turn, depends of the velocity of the CD and external density pex- 

Let Pi and P2 be velocities of the unshocked and shocked plasma in the frame of the 
forward shock. The shock jump conditions require [assuming pi = 0, ? ] 

^ (e2 - Pcx)(Pcx +P2) 
^2 ^ PJP^±P^ (40) 
62(62 -Pex) 

where €2 and p2 are energy density and pressure in the shocked medium and we assumed 
that the external medium is cold, pex = 0, eex = Pex- 

In the frame of the CD (CD is stationary with respect to downstream plasma), the shock 
has velocity V2, while the incoming flow has velocity vu- In the observer' s frame, the CD 
have velocity vcd, while the shock has 

J- + PCDP2 

Since the upstream plasma is at rest in the observer's frame 

PcD = V12 (42) 



Equations (42), an equation of state €2(^2), a condition of force balance on CD 

p, = = Ua{v = Vcd)'po/{2cx) (43) 
and shock jump conditions [9] 

/ gi+Pi+Pi V ^ ( S2 + P2+P2Y _ f £r+P.+P. ^ £2 + P2+P2 \ 

\ Pc. J \ P2 J \ PL P2 J 



constitute an equation on rjcD, the location of the contact discontinuity. (In Eq. (44) indices 



1 and 2 refer to unshocked and shocked external media correspondingly, pi = pex-) 
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For cold external plasma, pi = and assuming adiabatic index of F = 4/3 (so that 
^2 = 3p2 + P2), the shock jump condition gives 

P2 = ^ (jpi + 7pcx(48p2 + 49pex)) (45) 



Using continuity equations (40 42), solutions for the expansion flow given by (20) and 



shock jump condition (44) we find 



Po 



■>I3,CD) 



(46) 



Eq. (46) determines the Doppler factor of the contact discontinuity for plasma obeying 
equation of state with polytropic index F = 2, having initial Alfven Doppler factor 
expanding into a medium with density pex] the external medium is unmagnetized and obeys 



an equation of state with polytropic index F = 4/3. Eq. (46) is relativistically exact, no 
assumption about values of magnetization or external density were made, see Fig. |4j Note, 
that it correctly reproduces expansion into vacuum, pex = 0, 6yac = ^ao- Thus, shock jump 



conditions (44) are also applicable to the plasma- vacuum interface. 




Ci I Q. 




100 1000 



FIG. 4. Ratio of external density to the piston density Pex/po as a function of the Doppler factor 
of the contact discontinuity 6c d for one-dimensional expansion of strongly magnetized plasma 
with magnetization a, a = 100,10,1,0.1, top to bottom curves. Left Panel: Stationary initial 
conditions. Right panel: piston moving with 7^ = 10. Maximum values of 6cd = ^ao (^^^^ panel) 
and 6cD = ^ao^-w (right panel) correspond to zero external density, a vacuum interface. 



In the limit of relativistically strong shock, when the post-shock pressure is much larger 



than the post-shock density, relation (46) gives 

/3 PO 
= 32^^) 



3 B, 



2\ 1/4 



(47) 
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In this limit the velocity of the forward shock in the observer frame, Pfs = (Pcd + P2) /(I + 
PcdPi), is 

IFS = ^^cn =(l^A ={l^] (48) 



Pcx / \»P. 



cx , 



Heuristically, the expressions ( |47p8| can be obtained as a pressure balance at time t 



0+, when the contact instantaneously reach a Lorentz factor 7cd; the magnetic pressure 
in the plasma frame then is ~ 5^/7^^, and the kinetic pressure of the shocked external 
medium, ~ PexTcD- 

Note, that in the strong FS limit, the Lorentz factor of the CD, and of the forward 
shock, Eqns ( 47||48 ) are independent of cr, the composition of the driver. Only the value 



of total pressure in the piston matters (oc Bq in our case). Composition of the driver 
becomes important only for weak, non-relativistic FS, when the Lorentz factor of the CD 



(47) approaches the limit of expansion into vacuum. 7^ac = 1 + 2cr. This occurs for 



256 a 

In this limit of small external density, pex — 0, we find 



3 1 

Pex < TT^^Po (49) 



B. Expansion into medium: moving piston 



For moving piston, the system becomes somewhat more complicated. Similarly to the sta- 
tionary case, we expect the formation of four regions: undisturbed external plasma, shocked 
external plasma, undisturbed magnetized plasma; but the expanding magnetic piston may 
now turn into region of shocked piston materiel. These regions are separated by the forward 
shock, a contact discontinuity and an expansion wave or the reverse shock front. 

The cases, of the expansion wave and the reverse shock have somewhat different dynamics. 
In case of the expansion wave, the solution for simple waves derived in ^|TT] can be used to 
describe the expansion wave part of the flow. Alternatively, if RS shock is launched, which 
heats the piston material, the pressure on CD is determined by the sum of the magnetic 
pressure of magnetic field compressed at the RS and the kinetic pressure of shocked particles. 
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1. Forward shock and rarefaction wave 



In case of magnetic piston moving with initial Doppler factor 5^ into cold medium with 
density pex, the CD is moving with the Doppler factor (5^,cd,w determined from 



For strong forward shocks, the Lorentz factors of the CD and the FS are 

lCD,w '- 



(51) 





r 


, 8Pcx , 






1/4 


2pcx J 





= ( ] (52) 

(see Fig. |4]). 

Composition of the driver becomes important only for weak FS, when the Lorentz factor 



of the CD (47) approaches the limit of expansion into vacuum. F = 27^^,(1 + 2a). This 
occurs for 

The front of the rarefaction wave is stationary (located at ?7 = 0, 5,, = 1 when the flow moves 
sonically, with the Alfven velocity of the undisturbed plasma, (3^ = f3Afl- For higher the 
front of the RW is advected forward. 



2. Formation of a RS: supersonic motion of ejecta in the CD frame 



For sufficiently high velocity (3^, the location of the RW , Eq. (22), coincides with the 
location of the FS Sff^FS ~ '^Icd- This occurs for 



7„ > 27cDVCr 



(54) 



that is, when the wind velocity in the frame of the CD is supersonic, an obvious condition. 
In terms of initial magnetization and the ratio of densities, the RS forms when 



/3 Po _r2 



-a8 



A,0 



6^^3/2, a > 1 



(55) 



Relation (55) assume strong FS and high initial Lorentz factor 7^„ ^ 1. 
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3. Formation of a RS: subsonic motion of piston in the CD frame 



When motion of the piston in the frame of the CD is subsonic, yet the piston is moving 
faster than the CD, ■jcd < Iw < '^lcD\f^-, the flow is decelerated by a compression wave. 
For one- dimensional motion the compression waves are unstable to formation of shocks [? 



], so that the reverse shock will form for 7^ > 'Jcd, ond not at the condition (54)- As long 
as 7^ < 2'ycDy/o', the reverse shock is weak, in this range the RS may not form, if a more 
complicated flow patters are allowed. For •y^j > 2'ycDy/o', R-S becomes strong. (We define 
strong shocks as the shocks in which the upstream four-velocity in the frame of the shock is 
much larger than the upstream Alfven velocity). 



C. Reverse shock 

For sufficiently fast initial velocity, satisfying , 7^ 3> ■jw,crit, highly magnetized, a ^ 1, 
reverse shock forms. In this section we consider the dynamics of double-shock structures. 
For convenience we will make two approximations; first, we assume the the forward shock 
is strong and unmagnetized; second, we assume that the reverse shock is strong as well, 
Iw 3> ■yw,crit- We have to solve simultaneously for two shock jump condition and fiow 
continuity at the CD. 

In the frame of the CD, the FS is moving with P'pg = 1/3, while reverse shock is moving 
with [3] 

,,2 _ 1 + IQq + 8a^ + (1 + 2a)Vl + 16a + 16a^ 
Prs - iy ^ 26a + 8cr2 + (1 + 2a)Vl + IQcr + 16^2 ^ ^ 

In the frame of the shock, the post shock kinetic pressures are [3j 



plUf 4^272 



(57) 



while magnetic field satisfies 



5i7i/3i = (58) 



where indices 1 and 2 refer to quantities measured in the upstream and downstream in the 
frame of the shock. 

If the CD is moving with velocity (3cd and the initial velocity is the velocity of the 
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incoming plasma in the frame of the RS is 



/ _ (1 ~ I^cdI^w)I^'rs + I^w — (3cD /p-qN 

Using the above relations, the pressure behind the FS is 

(1 + 3fa)^ 4 2 

^^'^^ = 12(1 -fe) ^- " S^^^P- 

Kinetic pressure behind the RS is 

P2,,.n,RS = PO - g ^PO If ^ » 1 

(61) 



Magnetic pressure is 

P2,mag,RS = \iId11 i'l^CD {[^RS + /?») + PrsPu, + 1) ' ^ J^^O if > 1 (62) 



Balancing the pressures (60) and (62), we find 



-„„.(l,1.4.hj£)'^(^.jMy'%t.»l (63) 



The Lorentz factor of the CD in case of reverse shock coincides with the one obtained in case 



of rarefaction wave, Eq. (52). Thus, the forward shock is not influenced by the rarefaction 
wave-shock wave transition. We stress, again, that in the strong shock limit and for a ^ 1, 
the motion of the CD and the forward shock are independent of the composition of the 
piston (its density po)- 

If we express '-yen as a function of the total luminosity L/S = 7^(-Bo + Po) (assuming 

~ 1), where S is the cross-section of the flow, we find 

It is only weakly dependent on a. 

The RS is stationary in the observer frame when Pen = P'rs^ which for highly magnetized 
medium a ^ 1 (below, in this section, all relations are given for a ^ 1) gives 



7w = \L—(r (65) 
V 3 Po 

In this case 'Jcd = cr^^'^, which can be understood as the RS front is receding at almost 
with ua = o^l"^ in the CD frame. For higher the RS is advected towards the interface, 
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while for smaller 7^ it propagates in the opposite direction. Also, when --^cd = Iw, 1w 



(3/8)(po/Pex)o'5 there is no reverse shock or rarefaction wave; the motion of the wind 
matches exactly the motion of the CD. 



For values of 7^ sufficiently different from (65), the Lorentz factor of the reverse shock 
in the frame of stationary external medium is then 

1 / , 1cd\ I {usY^'^ ^ iftY^'^ ^c-^ ^ v^' moving forward 



iRS - \ - < 1/4 / \ 1/4" ^ 

2 V7CD / 1^ erfr-^ (^) if Icd < v^,RS moving backward 



(66) 



Recall, that relations (66) are applicable only if condition (55) is satisfied. 



D. Upshot: expansion of magnetized plasma into medium 

In Fig. |5] we qualitatively outline the dynamics of the fiow expansion in a medium. As 
one can see, the dynamics is rich in details and is much more complicated than for pure 
hydrodynamical expansion. Qualitatively, as the Lorentz factor of the piston increases for 
given density ration /, the rarefaction wave turns first into a weak RS and then into a strong 
RS. (In non-one-dimensional fiows weak reverse shocks can be avoided.) For extremely low 
external density, the expansion proceeds similarly to the case of expansion into vacuum. As 
the external density increases, first, the forward shock becomes strong. For higher p^x RS 
becomes backward propagating. 



IV. THREE DIMENSIONAL EXPANSION INTO VACUUM. 

In case of three-dimensional expansion of plasma carrying toroidal field, we expect that 
the flow will not be purely radial, as magnetic hoop stresses and magnetic pressure will, gen- 
erally, induce lateral motion (in 6 direction) in addition to radial expansion. Bur for strongly 
relativistic motion, corresponding to high magnetization parameter, the 6 component of the 
velocity will be much smaller, by a factor ~ 7^ than the radial component, so that the 
motion will be approximately conical, depending only on time and radial coordinate. 

For the three dimensional expansion, the governing equations are generically non-self- 
similar: there are two conserved quantities, energy and magnetic flux that scale differently 
with radius. Self-similarity can still be achieved in a narrow region, e.g. near the surface of 
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FIG. 5. Cartoon of the flow dynamics as a function of two parameters, 7^, the initial Lorentz factor 
for the piston and the ratio of piston density to external density / = High magnetization is 
assmned, a ^ \. Solid lines are related to the reverse shock, dashed lines are related to the 
forwards shock. The RS form for 7^ > 7cD) which is vTv^ \va.e. Below this line rarefaction wave 
propagates into piston. The RS shock is weak for jcD < < '^^CD\/o', the latter corresponds to 
"v/Jcr'^/^ line. In the frame of the undisturbed plasma the revers shock propagates in the forward 
direction for jcD > V^, which corresponds to -jw > \J f jo . The forward shock becomes non- 
relativistic for ^CD ~ 1) which is the Xj^fa] curve. The forward shock becomes weak (so that the 
post-shock temperature is non-relativistic) for ^CD ^ ^1 which is a 3/2/^ line. For weak FS, the 
expansion proceeds nearly as into vacuum. (Since relativistic motion is assumed, the lowest part 
of the plot, below the intersection of lines 7cz) = 1 and 7^ = 7cD, is outside of the applicability 
region.) 
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the bubble, where different radial scaling of energy and magnetic flux can be neglected. 

Similarly to the 1-D case, in case of three dimensional expansion of magnetized plasma 
into vacuum, the expansion front reaches terminal velocity immediately and coasts with 
constant velocity later. The constant terminal velocity in 3-D is the same as in 1-D - the 
vacuum interface always propagates with the terminal velocity 7 = 1 + 2(T, independent of 
the geometry [cf. [T8] . 

We expect that in the highly relativistic limit, the structure of the outflow will resemble 
a relativistic shock wave, where parameters change on a scale ~ 'j^^^ smaller than the overall 
size of the out flow. It is within this narrow region (and even smaller layer near the vacuum 



interface, see Eqns. (72 73])), that the self-similar solution derived below is applicable. ([18] 
did consider self-similar 3-D expansion for times much smaller than Tq/c, while astrophysical 
applications require times much larger than tq/c. ) 

For conical expansion of plasma carrying toroidal magnetic field, the conservation laws 
become 

9i(7p) + ^dr{r'j(ip) = 

dti^B) + -drinPB) = 
r 

dtToo + ^drir'Tor) = 

dtTor + ^drir^Trr) = (67) 



where we assumed a constant fixed polar angle 6 and neglected terms involving 6'-dependence. 

In three dimensions, the continuity and induction equations imply, B = [prBo) / (po'^o) 
where tq is an initial radius of the magnetized cavity. This scaling implies that for an outflow 
with a fixed energy, that starts with a finite region of non-zero magnetic field, one cannot 
assume constant magnetization, since a oc pr^. Then, the outflow dynamics is not self - 
similar, it depends on the initial o"(r). Self-similar solutions exist only very close to the edge 
of the expanding bubble; we derive them next. 

We assume that at the beginning of the expansion, close to the leading front, the plasma 
parameters (density and magnetic field) do not vary considerably. Eliminating density in 
favor of Alfven four- velocity p = [/^po/cr(ro/r)^, the condition that the magnetic and rest 
mass energy densities scale similarly with r requires Ua = g{x)/^- At the vacuum interface 
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the boundary condition is g = 0. Following [23], we introduce self-similar variable 

X={l + 2^lJ{l-r/t),x>l 



(68) 



and parametrize 7 = 7vac/- (Note that 7vac is a constant). In the leading orders in l/j^nc 
we find 

X//(1 + /)/' - gil - xfa + 2g'))g' = 

{l + xf)9f-2f{l-xf)g' = (69) 
The solutions satisfying boundary conditions /(I) = 1 and g{l) = are 



^1/3 



9 = (70) 

Energy densities per unit range of x are given by 

Usdx = (1 + n/2j'B' oc (1 - x-'^Ydx 

Updx = 7'P oc ix'^' - Iflxdx (71) 

These relations are applicable in the limit x ~^ 1- To estimate the region of applicability, 
we equate the total energy in the initial bubble, E ~ B^r^ to the energy within the narrow 
layer near the surface. We integrate Ub + Up from 77 = ! to 77 = ! + Ar] to determine Arj: 
E ~ B^rQi^Af])^ /ctq. Equating the two, we find Arj the width of the self-similar region in 
coordinate 77 

Ar/=(%o^)''' (72) 

In physical coordinates, 

where we used Fq ~ 2a. 

In conclusion of this chapter, we first reiterate that the terminal velocity of expansion 
into vacuum in 3-D case is the same as in ID: the terminal Lorentz factor is 7 = 1 + 2a 
(if starting at rest). The self-similar structure of the resulting bubble is applicable only in 
a narrow layer near the expanding vacuum interface, in most of the volume the field and 
velocity structure do depend on the particular initial conditions in the bubble. 
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V. DISCUSSION 



In this paper we found exact explicit solutions for one-dimensional relativistic expansion 
of polytropic fluid into vacuum and into plasma. In particular, we discussed an astropliysi- 
cal important case of strongly magnetized outflows; in this case especially simple analytical 
solutions can be obtained. We found exact solutions for one-dimensional expansion of mag- 
netized plasma into vacuum and into the cold medium both for stationary initial conditions 
and for a piston moving toward the interface. We found exact relations, applicable for ar- 
bitrary magnetization, relativistic motion and external densities. These results can be used 
for benchmark estimates of the overall dynamical behavior for the numerical simulations 
of relativistic plasmas, e.g. in heavy ion collisions, and in strongly magnetized outflows in 
particular. 

I am greatly thankful to Dimitros Gianios, Sergey Komisarov and Alexandre Tchekhovskoy. 
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